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We work out a semiclassical theory of shot noise in ballis- 



tic n + 



semiconductor structures aiming at studying 



two fundamental physical correlations coming from Pauli ex- 
clusion principle and long range Coulomb interaction. The 
theory provides a unifying scheme which, in addition to the 
current-voltage characteristics, describes the suppression of 
shot noise due to Pauli and Coulomb correlations in the whole 
range of system parameters and applied bias. The whole sce- 
nario is summarized by a phase diagram in the plane of two 
dimensionless variables related to the sample length and con- 
tact chemical potential. Here different regions of physical 
interest can be identified where only Coulomb or only Pauli 
correlations are active, or where both are present with differ- 
ent relevance. The predictions of the theory are proven to be 
fully corroborated by Monte Carlo simulations. 



I. INTRODUCTION 

Ballistic conductors are characterized by an active re- 
gion where carriers, once injected by contacts, move with- 
out suffering any scattering from contact to contact, i.e. 
the carrier mean free path is much longer than the sam- 
ple characteristic dimensions. In metals, ballistic trans- 
port is usually studied in point contacts. u Such con- 
tacts allow the analysis of elastic and inelastic scatter- 
ing processes-, by means of the so called point contact 
spectroscopyfl Since the Fermi wavelength in metals is 
very small (~ 0.5 nm), the nature of carrier transport 
is semiclassical and quantum effects related to the wave 
nature of the electrons can be disregarded. In semi- 
conductors, balii&tic transport has been investigated in 
bulk materials, pel point contactaJ and two dimensional 
electron gases;Q some studies have allowed the develop- 



ment of ballistic emission spectroscopy^ for the analysis 
of semiconductor heterointerfaces. Moreover, a variety 
of ballistic electrorL-cLevices with premising performances 
have been realizecffl or proposed. EjIIj Since the Fermi 
wave length in semiconductors can be as large as 40 nm, 
the nature of the carrier transport can be either semiclas- 
sical or quantum depending on the characteristic sample 
dimensions, L, with respect to the Fermi wave length, 
Xp. For Xp <C L transport is semiclassical, while for 
Xf ~ L transport is quantum. 

From a fundamental point of view, the study of non- 
equilibrium electronic noise (shot noise)EJ of ballistic con- 
ductors in the semiclassical transport regime offers a 
unique scenario where the simultaneous effect of two fun- 
damental physical interactions, namely the long range 
Coulomb interaction and the Pauli exclusion principle, 
can be investigated. On the one hand, Pauli correla- 
tions have an influence on the carrier injecting statis- 
tics at the contacts. On the other hand, Coulomb cor- 
relations can modify the carrier passing statistics inside 
the active regions. The effect of both mechanisms re- 
sults in values of the current noise below the full pois- 
sonian value, i.e. in shot-noise suppression. Thus, the 
non-equilibrium low frequency current spectral density 
is given by S/(0) = "f2ql, with I being the current flow- 
ing trough the sample and 7 a dimensionless parameter 
(Fano Factor) which takes values below unity because of 
the negative correlations induced by both mechanisms. 
The relative relevance of Coulomb and Pauli suppression 
in determining the Fano factor depends on parameters 
like the sample length, the temperature, the carrier den- 
sity or the applied voltage, thus allowing the monitoring 
of both correlation mechanisms. 

The first theoretical analyses on the non-equilibrium 
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noise properties of ballistic conductors were performed 
by Van der Ziel and .Bosnian in n + — rT — n + ballistic 
semiconductor diodese-a and by Kulik and Omel'yanchuk 
in ballistic metallic point contacts.E3 In the former case, 
by considering non-degenerate injection conditions, only 
the effect of the Coulomb correlations was evidenced. In 
particular, by using ideas borrowed from North's the- 
ory of vacuum tubes,E3 it was shown that under space 
charge conditions current noise results suppressed below 
the full Poissonian value (see also Refs. [ [l7],[l^]). In 
the later case, the high carrier density of metals pre- 
vented the presence of significant space charge effects, 
and hence the possibility to test the effect of the long 
range Coulomb correlations. In this case, by means of 
a selfconsistent semiclassical kinetic theory based on the 
Boltzmann-Langevin equation,E3 which included the ef- 
fects of the Pauli principle, it was shown that in the colli- 
sion free regime ballistic metallic point contacts displayed 
Nyquist noise with low frequency current spectral density 
given by Si(0) = AksT/Rs, where fcs is the Boltzmann 
constant, T the temperature and Rs = mvp/(q 2 n) the 
Sharvin contact resistance, with m being the electron ef- 
fective mass, up the Fermi velocity, q the unit charge, and 
n the carrier concentration inside the sample. This result 
is applied in the physical limit where qV, ksT <C Ep, and 
arbitrary relation between qV and ksT, where V is the 
applied bias and Ep the Fermi energy. These initial stud- 
ies did not allowed the study of the simultaneous effect 
of both Pauli and Coulomb correlation mechanisms. 

It was only very recently that ballistic structures in- 
cluding the simultaneous effect oL Pauli and Coulomb 
correlations have been analyzed£jt3 By making use 
of a semi-analytical approach,E3 Monte Carlo (MCI 
simulations ,Ej and an asymptotic analytical approach,E3 
non-equilibrium noise properties of ballistic semiconduc- 
tor structures have been investigated. It was shown that 
these structures, by allowing both space charge effects 
and degenerate injection statistics, were unique in or- 
der to study the simultaneous effect of both Pauli and 
Coulomb correlationsH3 However, the existing studies, 
being only valid under some limit conditions or per- 
formed for some specific sets of system parameters, do 
not offer a complete overview of the relative relevance of 
Pauli and Coulomb correlations in the whole range of sys- 
tem parameters and applied bias. This lack of a general 
overview can prevent the designing of suitable experi- 
mental strategies to test the theoretical predictions. The 
aim of the present work is precisely to address this issue. 

To this purpose we work out a general theory for 
the low frequency shot-noise properties of ballistic n + — 
i — n + semiconductor diodes that consistently describes 
also current-voltage characteristics. We consider ligthly 
doped bulk semiconductor materials since at low tem- 
peratures the bulk electrons are trapped by their par- 
ent donors, neutralizing them and thus minimizing both 
the e — e interaction and ionized impurity scattering, 
thus allowing long-jpacan free paths up to relatively high 
electron energies. ESEa The theory presented accounts for 



both the Pauli exclusion principle and the long range 
Coulomb interaction, and is applied to the whole range 
of system parameters (sample length, temperature, con- 
tact doping) and external bias. In particular, the theory 
allows us to study in a unifying framework the transition 
from: (i) nondegenerate to degenerate injection condi- 
tions, (ii) short to (asymptotically) long sample lengths 
and, (iii) low to (asymptotically) high applied voltages. 
By overcoming the limitations of the previous existing 
theories, the present study allows us to compress into a 
general scheme the full scenario displayed by the Pauli 
and Coulomb correlations in these ballistic structures. 
This scheme, validated and tested with a wide set of MC 
simulations, is believed to be of relevant assistance in de- 
signing future investigations on an experimental and/or 
simulation basis. 

The paper is organized as follows. In Sec. |l| we de- 
scribe the system under study. In Sec. |l| we present 
the physical model used to analyze the non-equilibrium 
noise properties. In Sec. IV the analytical expressions 



obtained for the transport and noise properties are re- 
ported, and their validity is checked by means of MC 
simulations. In Sec. ^ we propose a general scheme able 
to provide a systematic physical picture of the shot-noise 
properties of two-terminal ballistic conductors, and check 
its reliability by means of the developed theory. Finally, 
in Sec. ^ we summarize the main results of the paper. 
The appendix is devoted to technical derivations. 



II. SYSTEM UNDER STUDY 



The system under study is the n + 
semiconductor diode showed in Fig. |l|. 



n + ballistic 




x = 



x = L 



FIG. 1. Schematic band-diagram of the ballistic structure 
under study. 

It consists of an undoped semiconductor of length L, 
where injected carriers move in the total absence of scat- 
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tering events, sandwiched between two highly doped re- 
gions which act as ideal injecting contacts as specified 
later. The contacts and the ballistic region are taken 
of the same material (homodiode) or of different mate- 
rial (heterodiode). In the former case the band offset 
AE C ee 17(0) - E° c = U(L) - E£ (a symmetric structure 
is assumed for simplicity) vanishes, while in the later case 
it takes a finite value. The ballistic region is taken to be 
perfectly coupled with the contacts, thus no reflections 
take place at the interfaces. The contacts are assumed to 
be in quasiequilibrium at the given temperature, and at 
electrochemical potentials £o and with £t — £o = QV. 
The voltage drop in the contacts is assumed to be neg- 
ligibly small (Ohmic contacts), and hence all the band 
bending occurs in the active region of the sample. Under 
these assumptions, the effective contact chemical poten- 
tial fi, defined as fi = £o — U(0) = £l — U{L) is inde- 
pendent of the applied bias. Note that for the case of 
an homodiode fi coincides with the contact Fermi en- 
ergy, Ep = £o — E° = £o — while for the case of an 
heterodiode one has fi — Ep — AE C . We note the ex- 
istence of a maximum with amplitude U m in the poten- 
tial energy along the active region due to the presence 
of space charge. The structure is assumed to be suffi- 
ciently thick in the transversal directions so as to allow 
for a ID electrostatic treatment. Thus the-system is ID 
in real space and 3D in momentum space.113 The carriers 
are injected into the active region in accordance with the 
equilibrium Fermi-Dirac distribution function and sat- 
isfying the Pauli exclusion principle, thus .following the 
corresponding binomial injection statistics P 6 ^ 7 ! For sim- 
plicity, a single spherical parabolic band is assumed. 

The parameters describing the model system intro- 
duced above are therefore the length of the conductor, 
L, the cross sectional area A, the effective contact chem- 
ical potential, fi, the temperature, T, the static dielectric 
constant of the medium, e, the electron effective mass, 
m, and the charge carriers q. The independent-variable 
parameter is the applied voltage, V. 

As mentioned above, the correlations between differ- 
ent current pulses in the present system is twofold. On 
the one hand, Pauli correlations have an influence on the 
contact injecting statistics. Carriers above the Fermi en- 
ergy, at the tail of the Fermi distribution at the contacts, 
are injected obeying Poissonian statistics due to the neg- 
ligible influence of the Pauli principle at such high en- 
ergies. Carriers near and below the Fermi energy are 
injected following a binomial (sub-Poissonian) statistics 
due to the increasing influence of the Pauli principle at 
the highly occupied states. On the other hand, Coulomb 
correlations, through the potential energy maximum U m 
associated with the presence of space charge, can mod- 
ify the carrier passing statistics inside the active region. 
Indeed, the fluctuations of U m originated by the passage 
of carriers over the maximum modulate the transmission 
of further carriers (spacing them out more regularly in 
time) and smooth out the current fluctuations imposed 
by the random injection at the contacts. The relative 



importance of the two types of correlating mechanisms 
depends on the values of L, fi and V. 



III. PHYSICAL MODEL 

Within a semiclassical approach the description of 
the transport and noise properties of the system under 
study can be carried out by means of the Vlasov equa- 
tion, self-consistently coupled to the Poisson equation, 
and supplemented by appropriate fluctuating boundary 
conditions .113 In the one dimensional approximation fol- 
lowed here, the Vlasov-Poisson system of equations reads 
(in dimensionless units, as specified later): 



d d 8U(x,t) d 

\-v x — — 

dt dx dx dv x 



F(x,v x ,t) = 



d 2 U(x,t) 
dx 2 

n(x, t) 



n(x,t), 

SO 

dv x F(x 7 v x ,t) 



(1) 

(2) 
(3) 



where t is the time variable, x the spatial coordinate, v x 
the velocity in the x direction, U(x,t) the potential en- 
ergy, F(x, v x , t) the distribution function integrated over 
transversal momentum directions and n(x,t) the carrier 
density. In the Poisson equation the contribution of the 
intrinsic free carrier density is neglected in comparison to 
the injected carriers. Moreover, the boundary conditions 
for the potential energy read 



£7(0,i) = f7 = 0; U{L,t) = U L =V , 



(4) 



where we have assumed that the applied bias is fixed by a 
low impedance external circuit so that neither the electric 
potential nor the contact potential energies fluctuates, 
i.e., SU (t) = SU L (t) = 6V(t) = 0. Accordingly, the 
boundary conditions for the distribution function at the 
contacts read: 



F(0,v x ,t) = F(v x )+SF Q (v x ,t); 
F(L,v x ,t) = F(v x ) + 5F L (v x ,t); 



v x >0, (5) 
v x < . (6) 



These boundary conditions consist of two contributions. 
The first is of deterministic nature, and gives the aver- 
age value of the distribution function. The second is of 
stochastic nature, and describes the random injection of 
carriers. The average distribution function is given by: 



F{v x ) 



de±f FD (e± + e x - fi) 



ln(l 



) = fc (e x - fi) , 



(7) 



with e x — ^v x and e± = \v\ being the longitudinal and 
transversal kinetic energy, respectively, and where 
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Ifd 0) = 



1 



(8) 



is the Fermi-Dirac distribution function. Note that due 
to the integration over the transverse momentum direc- 
tions the effective one dimensional contact distribution 
function is f c (e) instead of Jfd (e)- Moreover, the fluc- 
tuating contributions SF (v x , t) and 8F^{v x , t) have. zero 
mean and low frequency spectral density given byEa 



5F a (v x ,t)5F a ,(v' x ,t')dt 
dfc (£x - M) 



d/j, 



-8 a , a 'S(e x -e' x )S(t-t') 



Ifd (e x ~ A 1 ) $a,a>8{e x - e' x )S{t - t') 



(9) 



Equations (0)-(||) constitute the complete set of equa- 
tions to study the noise properties of the ballistic struc- 
ture described in Sec. |l] 

In the previous set of equations, and in what follows, 
we use dimensionless variables to simplify the notation. 
The dimensionalizing factors used for energy, length, car- 
rier density, velocity, distribution function, electric po- 
tential, electric field, electric current, resistance and cur- 
rent spectral density are, respectively, 



TT J, T J €kBT N A (*T) 3 / 2 



q 2 N 



h 3 



(k B T\ 1/2 N k B T 

v = ;F = — ;V = 

\ m J v q 



(10) 



E = ^ N ° k ^ T ; l Q = qAN v ; R = ; Sj = 2ql , 

where h is the Planck constant. With these dimension- 
alizing factors the dimensionless parameters describing 
the system are L/Lq — > L and [i/k B T — » /i, and the di- 
mensionless independent variable is qV/k B T — > V. The 
physical meaning of the different parameters is the follow- 
ing: iVo is the effective density of states in the conduction 
band, L is the Debye screening length associated to No, 
Uq is the thermal energy, w is the thermal velocity, Vo is 
the thermal voltage, Eq the thermal electric field, Jo the 
current associated to iVo and vq, Rq the resistance asso- 
ciated to Vq /Iq and 5/ the shot noise level associated to 
A- 

To complete our analysis, we will perform MC simula- 
tions of the system under study, so that the validity of 
the analytical theory can be checked by comparison. To 
this end we use an ensemble MC simulator, 3D in mo- 
mentum space and self-consistently coupled with a ID 
Poisson solver to account for Coulomb interaction. The 
carrier dynamics is simulated in the ballistic active re- 
gion of the structure and the electron injection from the 
thermal reservoirs is modeled according to Fermi statis- 
tics. Due to Pauli principle, the instantaneous occupancy 
of an incoming electron state with energy e and imping- 
ing at the interface between the ideal thermal reservoir 



and the active region fluctuates in time obeying a bino- 
mial distributiorc3 with a probability of success given by 
}fd{s — This statistics is implemented in the MC 
simulation of the contact injection by introducing a dis- 
cretization of momentum space and using the rejection 
technique to select the times of injection at every mo- 
mentum state.EH As limiting cases, when e — /i <C — 1, 
/fd{£ — A 4 ) — 1 and the injection statistics of the cor- 
responding state is uniform in time. By contrast, when 
e — fi » 1, /j?b(£ — ju) <C 1 and the injection statistics is 
Poissonian. For the calculations we use the following pa- 
rameters: T — 300 K, m = 0.25 uiq, e = 11.7 eo, with uiq 
the free electron mass and eo the vacuum permittivity. 



IV. ANALYTICAL SOLUTION 

The low frequency solution of the model presented in 
Sec. [II can be obtained in fully analytical form. For 
the sake of conciseness, below we only present the final 
expression of the relevant quantities of interest and refer 
to the appendix |X| for the details of the derivation. 



A. I — V characteristics and steady-state profiles 

Following the results presented in appendix |a|, the cur- 
rent voltage (I — V) characteristics of the ballistic con- 
ductor presented in Sec. || can be calculated as 



du [f c (u - fo) - f c (u - fo)] 



(11) 



where / c (e) is given in Eq. U m is the average value 
of the maximum potential energy at applied bias V, and 
£i = V + /x and £o = M are t ne electrochemical potentials 
at the contacts located at x — L and x — 0, respectively. 
The value of U m can be calculated through the following 
equation (see appendix [A|) 



L = 



dU 



dU 



u E-(U,U m ) Ju L E+(U,U m ) 
where 

E-(U,U m ) = 



(12) 



du 



U — U r , 



+ J duVu - U2f c {u - Co) 
and 

E+(U,U m ) = 



[/c(u-&) + /c(u-&)] 

2* (13) 



du 



y/u — U — \ U — U 



[/c(«-&)+/c(u-&)] 
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+ / duVu zr U2f c (u - £ L ) } 23 




Figure || shows the I — V characteristics as calcu- 
(14) lated from Eq. ( [TT| ) for fx — 8 and several sample lengths. 
The figure also reports the results of MC simulations. 
We have found perfect agreement between the results 
of the theory and those of MC simulations as expected, 
since both results are exact solutions of the same physi- 
cal model. The main features of the I —V characteristics 
consist in the presence of a linear (Ohmic) behavior at 
low applied voltages and a current saturation regime at 
high voltages. The resistance of the linear regime is given 
by (in dimensionless units) 
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FIG. 2. Maximum energy potential Um/ksT as a function 
of voltage for jit/fcsT = 8 and several values of the sample 
length L/Lq. Solid line refers to the result of the analytical 
calculations and open circles to those of the MC simulations. 

Figure |^ reports the maximum potential energy as a 
function of the voltage, as calculated from Eq. (|12|) with 
/i = 8 at different values of the sample length. MC re- 
sults are also shown. An excellent agreement is found be- 
tween analytical and MC calculations. At vanishing volt- 
ages U m saturates at higher values the longer is the nor- 
malized length, because of the increasing effect of space 
charge inside the ballistic region. At increasing voltages 
U, n starts decreasing till vanishing at a threshold voltage 
V m . 
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qV/k B T 

FIG. 3. I-V characteristics for \jLjksT = 8 and several val- 
ues of the sample length L/Lo. Solid lines refer to the re- 
sult of the analytical calculations and open circles to those 
of MC simulations. Dashed curves indicate the typical power 
law behaviors of the Ohmic (V) and Child-Langmuir (V 3 ^ 2 ) 
regimes. 



(15) 



In general, this value depends on both the sample length 
L and the reduced chemical potential ll. When U m <C \li\ 
the resistance takes a value independent of the sample 
length K %™ v = l/f c (—fi) in agreement with the results 
of Ref. [ E8fl. Typically this happens for L <C 5 and 
/i> 5, which represent the conditions for negligible space 
charge inside the ballistic region. In this limit, and when 
fi ^> 1, Re™ zv tends to the Sharvin contact resistance 
Rs- The value of the saturation current is given by (in 
dimensionless units): 



f c (u - fi)du . 



(16) 



Current saturation takes place when all carriers injected 
from one contact reach the opposite one, while none of 
the carriers injected from the other contact is able to 
cross the ballistic region. The value of the saturation 
current is independent of the sample length since it is 
only determined by the emission properties of the con- 
tacts. The voltage for the onset of current saturation, 
V cr , depends on both L and fx, as can be seen in Fig. ||. 
In the case that // > V™, where V™ is the bias value 
at which the maximum of the potential energy disap- 
pears (see Fig. |J), one has V cr = fx + 3. Otherwise, when 
ll < V™, one has V cr — V™, and saturation coincides 
with the disappearance of the maximum of the poten- 
tial energy. Between the linear and the current satura- 
tion regimes the I — V characteristic displays a nonlin- 
ear region, whose properties are determined by the pre- 
cise values of L and fj,. The asymptotic behavior of this 
non-linear region for large values of L in the limit when 
U m <C V < V cr consists of a slightly sublinear region fol- 
lowed-iby a superlinear characteristic of Child-Langmuir 
typeEa I CL = V 3 / 2 /L 2 , as has been found in Ref. [ p2[ . 
The sublinear region is associated with a non zero tem- 
perature value which smooths out the Fermi distribution. 
The superlinear region is associated with space-charge ef- 
fects driven by the presence of U m . 

The steady state profiles can also be calculated in 
closed analytical form. The potential energy profile can 
be obtained in inverse form from the following relations 
(see appendix |A|) : 
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,U (x) 



dU 



Uo E- (U, U m ) 
dU 



< x < x r , 



(17) 



U + (x) 



Ul E+(U,U m ) 



= x — L 



x m < x < L , (18) 



where x m is the location of the potential energy max- 
imum inside the ballistic region. The value of Xfn Call 
be calculated from either Eq. ( p"7j ) or Eq. ( |l8| ) by setting 
the value of U(x) equal to U m . Once the potential en- 
ergy profile is obtained, the electric field profile can be 
calculated as (see appendix |A]) 



E(x) 



E-lU (x),U, 
E+ (C7 + (ar),C7„ 



< x < x m 
x m < x < L 



, (19) 



where E~ (U, U m ) and E + (U, U m ) are given in Eqs. 
and ( |l4|) , respectively. Finally, the carrier density profile 
is obtained from (see appendix ^) 



n(x) = 

( r+oo 
J~U m 



JU (x) 



du 



2(u-U (xy 



du 



J U + (x) 



2(u-U (x) 
du 



[/c(u-£o) + /c(u-£ L )] + 
= 2/ c (u-^o); 0<x<x m 



2(u-U (x)) 



(20) 



da 



2(u-U + (x)) 



=2/ c (u-^i); x m < x < L 



Figure [l] reports the steady state profiles for \i = 8 
and L = 8, at several values of the applied bias. For 
the sake of comparison, the same figure also reports the 
profiles obtained from MC simulations. Again analytical 
calculations are found to perfectly agree with the sim- 
ulations. The figure clearly illustrates the space-charge 
nature of the transport and the presence of the poten- 
tial energy maximum. The value of the potential energy 
maximum decreases systematically at increasing bias (as 
already indicated in Fig. while its location shifts to- 
wards the left contact. At the same time, the free charge 
redistributes inside the structure by shifting the mini- 
mum of the carrier density from the center of the sample 
towards the right contact. 



B. Noise properties 

The low frequency noise properties of the system un- 
der study are characterized by the low frequency spectral 
density of current fluctuations defined as (in dimension- 
less units), 



Si(0) 



6I(0)SI(t)dt 



(21) 



When interested in the study of the microscopic corre- 
lations the proper figure of merit is the Fano factor 7, 
which is obtained from 5/(0) as (in the dimensionless 
units) 
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FIG. 4. Steady state profiles of the potential energy 
U(x) jksT ', electric field E(x)/Eq and carrier density n(x)/No 
for n/kBT — 8 and L/Lq = 8, at different values of the ap- 
plied voltage qV/ksT. Solid lines refer to the results of ana- 
lytical calculations, open circles to those of MC simulations. 

For conditions far from thermal equilibrium (V > 3) a 
value of 7 = 1 corresponds to the absence of correlations 
between different current pulses, while 7 < 1 (7 > 1) 
corresponds to the presence of negative (positive) corre- 
lations. 
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As shown in appendix |X[ 5/(0) can be calculated in 
closed analytical form for the model presented in Sec. HI. 



The final result can be written in the following compact 
form (see appendix ^) 



r+oo 

5/(0) = /_ du [ lL {u)\ 2 f FD (u - + 
Ju L 



(23) 



du[j (u)] /fd(u-Co) 



u 



where (see appendix |A|) 



1l{u) 
7o (u) 



U m < u < +oo 
U L < u <U m 



1 + Qj>(u) 
fi7<(tt) 

-l + fi ; 7 > (u) Z7 TO < u < +oo 
U < u < U m 



Here, we have defined 



Q = — 



(24) 
(25) 

(26) 



with 



Eq El 



dwy£{u)f FD (u-£ L ) 



+ 00 



Moreover 



+ /_ du^>{u) [f FD (u - Co) + f FD (u - &)] • (27) 



U m ^2 (u - U) - J 2 (u-U m ) 
dU— V 



E- {U,U m f 



u m y/2 {u - U) -\/2(u- U m ) 
dU 



Ur 



E+ (U,U m ) 3 



t£(«) 



-2/ W^Z?\ 
lu L E+ (U, U m ) 



7o<(u)=2 / dU 



u E~(U,U n 



(28) 
(29) 
(30) 



From the previous equations we can evaluate 57(0) and, 
in turn, the Fano factor 7. Note that these analytical 
expressions are valid in the whole range of system pa- 
rameters L and [i, and in the whole range of applied bias 
V. 

For the purpose of a reliability test, the results ob- 
tained from the analytical formulae are compared with 
those of MC simulations in Fig. |j| Here, the low fre- 
quency current spectral density is reported as a function 
of applied bias for \i — 8 and several sample lengths. 
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FIG. 5. Low frequency current spectral density Si(0)/2qIo 
as a function of the applied voltage qV/ksT for /i/fcsT = 8 
and several values of the sample length L/Lq. Solid lines refer 
to the results of analytical calculations, open circles to those 
of MC simulations^ Dot-dashed lines are the results of the 
asymptotic theoryE3 

As can be seen in the figure, the agreement between 
the analytical theory and MC simulations is excellent, 
thus proving the reliability of both theory and simula- 
tions. For the sake of completeness Fig. || also reports 
the results obtained by means of the asymptotic theory 
developed in Ref. [ ^2j. According to this theory, in the 
limit of U m <C V < V cr the low frequency current spec- 
tral density can be approximated by (in dimcnsionlcss 
variables) 



with 



Sj sym = p(p - U m ) =- 

1 ^ J v + u„ 



K ' \ 4F (o)Fi(a) 



(31) 



(32) 



where Fj(a) = 1/T(j + 1) dyy 1 f FD {y - a), with 
r(z) being the Gamma function. As seen in Fig. I the 
asymptotic theory agrees well with the present theory 
for U m <C V < V cr - However, it cannot describe ei- 
ther the transition between low and high bias regions, or 
the case of small sample lengths in which the condition 
Um <§C V < V cr can not be satisfied. 

The main features of the low frequency spectral den- 
sity depicted in Fig. [| are summarized as follows. At 
low voltages (V < 3), where Ohmic conditions are sat- 
isfied, 5*7(0) is bias independent and takes the value (in 
dimcnsionless variables) 



5f (0) = — 



(33) 



in agreement with Nyquist theoremS At high voltages 
(V > V cr ), when the sample exhibits current saturation 
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conditions, £7(0) is again bias independent and given by 
(in dimensionless variables) 



Sf at (0) = 



+ OC 



duf FD {u - £ ) = /c(~m) 



(34) 



In the region of intermediate values of voltages, <S7(0) 
can present a monotonic or a non-monotonic behavior 
(with the presence of a minimum) determined by the in- 
terplay between Coulomb and Pauli correlations, as will 
be detailed in the next section. 



V. THE SCENARIO OF SHOT-NOISE IN 
BALLISTIC N + — J — N + STRUCTURES 

The present theory enables us to investigate separately 
the relevance of the two mechanisms responsible for the 
correlations in the system under study, namely, the Pauli 
exclusion principle and the long range Coulomb interac- 
tion, in the whole range of system parameters. Accord- 
ingly, we propose a general scheme which summarizes the 
whole scenario of the shot-noise properties exhibited by 
n + — i — n + ballistic semiconductor structures. To con- 
struct such a scheme, we make use of the Fano factor, 
7, which is factorized into the two independent contri- 
butions, 7p and 7c, related to the Pauli and Coulomb 
correlations, respectively. Indeed, according to Eq. (||), 
the fluctuations of the contact distribution function at 
different energy levels are uncorrelated, so that the only 
source of correlations among carriers injected with dif- 
ferent energy is the Coulomb interaction in the active re- 
gion. As a consequence, both contributions to the Fano 
factor are independent, which implies 7 = "fpjc- Thus, 
the Pauli contribution jp corresponds to the Fano fac- 
tor that would be obtained in the absence of the self- 
consistent long range Coulomb interaction. It can be 
easily evaluated from the noise calculation performed in 
the appendix [a] by neglecting the self-consistent contri- 
bution. This is equivalent to set up in Eq. (|23 



7o (u) 



1 U m < u < +00 

U L < U < U m 

— 1 U m < u < +00 
U < u < U m 



(35) 
(36) 



For the current spectral density associated with Pauli 
correlations only, it is thus obtained 



Sf(0) = f e (U m - &) + f c {U m - Co) 



(37) 



so that the Pauli contribution to the Fano factor is found 
to be 



Sf(0) = MU m -tL) + MU m -Zo) 
I " g°°du[f c (u- 



IP 



where use is made of Eq. (O 



(38) 



The Coulomb contribution is then evaluated as 7c = 
7/7P. According to these definitions, for V > 3 values of 
7p 7^ 1 correspond to the presence of Pauli correlations, 
while values of 7c ^ 1 correspond to the presence of 
Coulomb correlations. 

A detailed analysis of Eq. (pq) indicates that the Pauli 
contribution is essentially dependent on the difference 
U m — [x. Hence, when U m — fx > it is 7p — > 1, thus indi- 
cating the absence of Pauli correlations. In particular, for 
nondegenerate injection conditions jx < 0, one always has 
that U m — fx > 0, hence indicating the absence of Pauli 
correlations, as should be. By contrast, when U m — Li < 
it is 7p < 1, thus indicating the presence of Pauli cor- 
relations. Since U m is a decreasing function of the bias, 
the condition U m — fx < implies automatically that the 
inequality U m — ll < is satisfied for all bias values, and 
hence the presence of Pauli correlations for V > 3. On 
the contrary, when U m — /1 > Pauli correlations are 
absent for low (or intermediate) bias values and present 
for bias values sufficiently high to validate the condition 
U m - n<0. 

Concerning Coulomb correlations, their presence or ab- 
sence is roughly determined by a value of the ratio L/Ld c 
higher or lower than unity, respectively, where L d c is the 
Debye screening length corresponding to an homogeneous 
system with charge density equal to the equilibrium con- 
tact density. Taking into account the effects of degener- 
ancy, Lp> a is calculated as (in dimensionless variables) 



L, 



( dn e c 



-1/2 



1 



2 J +o ° du 



frp("-c) 

2u 



(39) 



where we used that from Eq. (pOj) one has 



'(0) 



+ °° a tab 

du -= 

11 \j2u 



(40) 



On the basis of these considerations, for the scenario 
of the shot-noise properties in ballistic conductors we 
propose the general scheme displayed in Fig. ^. In this 
scheme we identify five different regions in the plane (L, 
Li), corresponding to five different possibility of interplay 
between Coulomb and Pauli correlations. The different 
regions in Fig. [(] are determined by the three lines defined 
by the equalities L = Ld c , (J, = and ll = U^. 

The shot noise behavior in each of the five regions cor- 
responds, respectively, to: (i) the absence of both Pauli 
and Coulomb correlations {ll < and L/Lp> c < 1); (ii) 
the presence of Pauli correlations and the absence of 
Coulomb correlations {ll > 0, — ll < and L/Ld c < 
1); (iii) the absence of Pauli correlations and the presence 
of Coulomb correlations {fj, < and L/Lp> c > 1); (iv) the 
absence (presence) of Pauli correlations for low (high) 
bias and the presence of Coulomb correlations (/1 > 0, 
U m — /i > and L/Lp> c > 1); (v) the presence of both 
Pauli and Coulomb correlations {fx > 0, U m — /1 < and 
L/L De > 1). 
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FIG. 6. Parameter plane in the length /chemical poten- 
tial space representing the five different behaviours of the 
shot-noise properties of a ballistic conductor according to the 
relevance of the Pauli and Coulomb correlations. The three 
continuous lines define, respectively, the equalities: L — Ld c > 
H = and = fi. 

The reliability of the above scheme has been tested by 
performing a series of theoretical calculations for the rel- 
evant regions identified above. We have found that the 
proposed scheme is essentially valid, except in the zones 
close to the lines separating the different regions, where 
intermediate behaviors have been observed. Represen- 
tative examples concerning the Fano factor, 7, and the 
contributions into which it is decomposed, jp and 7c, 
are shown in Figs. |?] to |l^ for each of the five regions 
individuated in the L-[i plane. 



10 1 



10° 




u/k B T = -5, L/L = 5 
region (i) 



10 



u/k B T = 8, L/L = 0.5 
region (ii) 



10 1 



10° 1 1 

qV/k B T 



10 2 



FIG. 7. Fano factor as a function of the applied volt- 
age qV/ksT for L/Lo = 5 and /i/fegT = —5 correspond- 
ing to region (i) of Fig. ^| (solid line) and L/Lo = 0.5 and 
/j/feT = 8 corresponding to region (ii) of Fig. || (dot-dashed 
line). Dashed line represents the Coulomb contribution to the 
Fano factor. By definition, in this case the Pauli contribution 
to the Fano factor is indistinguishable from the actual Fano 
factor. 



Figure ^ displays the Fano factor for L = 5 and \i = — 5 
(L/Ljj o = 0.32) corresponding to region (i) (continu- 
ous line) and for L = 0.5 and /i = 8 (L/Lp, c = 0.71) 
corresponding to region (ii) (dot-dashed line). In both 
cases the presence of space charge in the active region 
is nearly negligible, insufficient to originate a potential 
energy maximum large enough so as to lead to Coulomb 
suppression. Therefore, both regions concern with the 
absence of Coulomb correlations (jc = !)• In region (i) 
Pauli correlations are absent (<yp = 1 for V > 3), since 
fi = —5 implies that carriers are injected at the contacts 
with energies e — fi > 3, so that the distribution function 
at the contacts is well approximated by the nondegen- 
erate Maxwell-Boltzmann distribution, and the injection 
statistics is poissonian. On the contrary, in region (ii) 
Pauli correlations are responsible for the suppression of 
shot noise (jp < 1, for V > 3) since degenerate injection 
conditions prevail. 



10 1 



10° 



10-" 







u/k B T = -3, L/L = 40 - 






region (iii) 






■ . Y P 






\\ / '. 

V / '. 



10" 



10° 



10 1 



10 2 



qV/k B T 



FIG. 8. Fano factor as a function of the applied voltage 
qV/ksT for L/Lo = 40 and fi/ksT = —3 corresponding to 
region (iii) of Fig. ^| (solid line). Dashed and dotted lines 
represent, respectively, the Coulomb and Pauli contributions 
to the Fano factor. Dashed-dotted line represents the results 
of the asymptotic theory. 

It is worth noting that, in all cases where the long 
range Coulomb correlations are absent, as in the present 
case, the Fano factor at low voltages decreases inversely 
with the applied bias according to the law 



1 V ' 



(41) 



which corresponds to the thermal noise behavior, while 
at high voltages (V > V cr ) it becomes constant with a 
value given by 



Zc(-m) 



Jo °° fc(u - ^) du 



(42) 



which corresponds to current saturation conditions. To 
obtain Eq. (E|) use is made of Eqs. (|l6|) and (EmI). The 







value of r ) sat interpolates monotonically between 1 for 
fi < —3 (nondegenerate injection statistics) and 2/ji for 
fj, > 3 (strongly degenerate injection statistics). Since 
the Coulomb correlations vanish identically at thermal 
equilibrium and under current saturation conditions, the 
two limiting behaviors represented by Eqs. (|4l|) and 
are common to all cases, as we will see below. 

Figure H reports the Fano factor for L — 40 and fi = — 3 
(L/Ld c =6.9) corresponding to region (hi) of the general 
scheme. It shows the presence of Coulomb correlations 
(7c < 1 for 3 < V < V cr ) and the absence of Pauli 
correlations (jp = 1 for V > 3). Pauli correlations are 
absent because \i = — 3 < (see above), while Coulomb 
correlations set in due to relevant space charge effects. 
Due tq-fche small value of the ratio L /L d c the asymptotic 
theoryca provides only some rough agreement with the 
exact result in the limited range of voltages associated 
with a suppressed behavior. For larger values of the ratio 
L/Ljj c the agreement improves. The behavior displayed 
by the system in region (hi) was previously examined in 
detail in Ref. [£§]. 

Figure |^ reports the Fano factor for L — 20 and ji = 8 
(L/Ld c = 28.2) corresponding to the region (iv) of the 
general scheme. Here we assist to a low and intermediate 
voltage range where Coulomb correlations are dominat- 
ing and to a high voltage range where Pauli correlations 
prevail. Indeed, the Fano factor displays the presence of 
Coulomb correlations (7c < 1 for 3 < V < V cr ), the ab- 
sence of Pauli correlations in an intermediate bias range 
(7p ~ 1 for 3 < V < 10), and the presence of Pauli 
correlations for higher bias (■yp < 1 for V > 10). 
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FIG. 9. Fano factor as a function of the applied voltage 
qV/ksT for L/Lo = 20 and /i/fc_gT = 8 corresponding to 
region (iv) of Fig. ^ (solid line). Dashed and dotted lines 
represent, respectively, the Coulomb and Pauli contributions 
to the Fano factor. Dashed-dotted line represents the results 
of the asymptotic theory. 

Physically, this behavior is due to the fact that in this 
region, even if contacts are under degenerate injection 
conditions (fi > 0), one has f/^, — /1 > 0, so that for the 



lowest applied voltages (3 < V ^ 10) it is U m — \i > 
and the current flows is only due to Poissonian carri- 
ers at the tail of the Fermi distribution function, lead- 
ing to 7p ~ 1. As V increases, it is U m — fi < and 
sub-Poissonian carriers near and below the Fermi level 
increasingly contribute to the current and low-frequency 
noise, so that Pauli correlations become manifest in the 
noise, thus leading to jp < 1. Note that in region (iv) 
one has generally JJLp> a ^ 10, thus ensuring that the 
asymptotic theorjo provides a good approximation to 
the exact theory presented here in the range of bias sat- 
isfying U m <C V < V cr , as illustrated in Fig. |[ 

Finally, Fig. |l^ reports the Fano factor for L = 2 and 
(X = 8 (L/Lp> c = 2.8) corresponding to region (v) of the 
general scheme. The Fano factor displays the presence of 
both Pauli correlations (jp < 1 for V > 3) and Coulomb 
correlations (7c < 1 for 3 < V < V cr ). In this region 
one always has U rn — fi < thus ensuring that Pauli 
correlations are present at all applied bias. Therefore, 
this is the most interesting region to analyze the interplay 
between Coulomb and Pauli correlations. We note that 
here the exact theory presented in this paper becomes 
strictly necessary, since in this region the ratio L/Lp> c 
never takes values much higher than unity (one typically 
obtains 1 < L/Ld c ^ 10), thus reducing significantly the 
usefulness of the asymptotic theoryj23 as illustrated in 
Fig. 0. 
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FIG. 10. Fano factor as a function of the applied voltage 
qV/ksT for L/Ln = 2 and fi/ksT = 8 corresponding to re- 
gion (v) of Fig. (solid line). Dashed and dotted lines rep- 
resent, respectively, the Coulomb and Pauli contributions to 
the Fano factor. Dashed-dotted line represents the results of 
the asymptotic theory. 

To provide more insight into the interplay between 
Coulomb and Pauli correlations in region (v), Fig. |ll| 
reports the Fano factor as a function of bias for L = 3 
and several values of fi belonging to region (v). It is ob- 
served that by increasing the value of \x the contribution 
of the Pauli correlations decreases faster than that of the 
Coulomb correlations. This is understood by noting that 
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for /i > 3, 7p varies as 2//x, while 7c, being determined 
by the ratio L/Lo a , varies as y/JI. Physically, this be- 
havior reflects the increasing amount of sub-Poissonian 
carriers contributing to the current as \i becomes higher. 
These carriers lead to a more and more pronounced Pauli 
suppression while they do not contribute so much to 
Coulomb suppression, since due to their sub-Poissonian 
character the fluctuations of U m they induce are much 
less significant than those originated by the Poissonian 
carriers above the Fermi level. For this reason, jc tends 
to saturate for the highest values of V < V cr while 7 P 
continues decreasing. Remarkably enough, in the region 
(v) for V > 30 the total Fano factor 7 is well described 
in all cases by a 1/V dependence. As seen in the fig- 
ure, this dependence on voltage is due to the joint action 
of both Pauli and Coulomb correlations. It can also be 
observed that the onset of Coulomb suppression takes 
place for higher values of V as [i increases, since for this 
suppresion to become significant it is necessary that the 
contribution to the current due to carriers injected at the 
right contact becomes negligible. This takes place when 
V + U m = fx, condition that requires a higher value of V 
as fi increases, for a given value of L. 
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FIG. 11. Fano factor as a function of the applied voltage 
qV/ksT for L/Lq = 5 and several values of fi/ksT corre- 
sponding to region (v) of Fig. |^ (solid line). Dashed and dot- 
ted lines represent, respectively, the Coulomb and Pauli con- 
tributions to the Fano factor. The dashed-dotted line gives 
the 1/V slope. 

To illustrate a realistic example where the theoretical 
results presented here could be applied, we consider as a 
significant example the case of a ballistic GaAs n + —i—n + 
homodiode at T — 4 K and contact density n c = 1.14 
10 16 cm" 3 . For this structure k B T/q = 0.3 mV, fi = 8 
(0.24 meV), L ~ 31 nm and L Dr ~ 0.7 (22 nm), where 
we have taken m = 0.066 mo, and e = 12.9 eo. For a 
sample length of 500 nm, L ~ 16, and according to Fig. ^ 
we are in region (iv) of the general scheme. For a sample 
length of 70 nm, L ~ 2.3 and we are in region (v). For 
a sample length of 15 nm, L ~ 0.5 and we are in region 



(ii). To explore regions (i) and (iii) one can refer to the 
case of heterodiodesH We conclude that in principle the 
different behaviors of the Fano factor predicted by the 
present theory can be investigated experimentally within 
realistic conditions. 



VI. CONCLUSIONS 

We have presented a semiclassical theory of non- 
equilibrium noise (shot-noise) properties of n + — i — n + 
ballistic semiconductor structures aiming at evaluating 
the relevance of the long range Coulomb interaction and 
of the Pauli exclusion principle. The theory covers the 
whole range of system parameters, physically identified 
in the contact chemical potential, the sample length, and 
the applied voltage. Within a unitary scheme free from 
any approximation we succeed in investigating the tran- 
sition between: nondegenerate and degenerate injection 
conditions, short and long sample lengths, and low and 
high applied bias. Through the determination of the 
Fano factor, we have analyzed the relevance of Coulomb 
and Pauli correlations. At applied voltages above about 
ZksT/q both correlations lead to the suppression of shot 
noise. We have identified five different regions in the 
plane defined by the sample length and the chemical po- 
tential corresponding, respectively, to the following con- 
ditions, (i) the absence of Pauli and Coulomb correla- 
tions; (ii) the presence of Pauli correlations and the ab- 
sence of Coulomb correlations; (iii) the absence of Pauli 
correlations and the presence of Coulomb correlations; 

(iv) the absence (presence) of Pauli correlations for low 
(high) bias and the presence of Coulomb correlations; and 

(v) the presence of both Pauli and Coulomb correlations. 
Case (i) corresponds to small sample lengths and non- 
degenerate conditions when different current pulses are 
clearly uncorrelated. Case (ii) occurs for small sample 
lengths and degenerate injection conditions so that Pauli 
correlations are the only being active. Case (iii) implies 
a large sample length and non-degenerate injection con- 
ditions so that Coulomb correlations are the only being 
active. Case (iv) is associated with large samples and de- 
generate injection conditions. Here Coulomb correlations 
are always present because of the small value of the Debye 
screening length. In addition, at low bias only Poissonian 
carriers at the tail of the contact energy distribution con- 
tribute to the noise so that Pauli correlations are absent, 
while at higher bias Pauli principle becomes active due 
to the increasing amount of carriers obeying binomial in- 
jection statistics that contribute to noise. Finally, case 
(v) refers to moderately long samples and degenerate in- 
jection conditions so that both Coulomb and Pauli cor- 
relations are present simultaneously. The results of the 
theory are in perfect agreement with analogous MC sim- 
ulations. Therefore, besides offering a complete physical 
picture of the subject, this study provides new insight 
into the noise properties of ballistic conductors, and of 



11 



mesoscopic systems in general. We believe that the the- 
ory here developed constitutes a powerful tool to design 
experimental investigations of the nonequilibrium noise 
properties of solid-state ballistic conductors. 
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APPENDIX A: DERIVATION OF THE MAIN 
FORMULAE 

In this appendix we derive the main formulae used in 
the paper by following a procedure similar to that used 
in Ref. [ [l8| , here extended to the case of degenerate in- 
jection conditions. To this end, the first step consists 
in solving for the distribution function that satifies the 
Vlasov-Poisson-Langevin system of equations, Eqs. (|l|)- 
(ph. To find the solution we first note that, since trans- 
port is ballistic, the total longitudinal energy of a carrier 
u, defined as 



1 2 

2 * 



U(x,t) 



+ U(x,t) 



(Al) 



remains constant during the flight of the carrier through 
the structure. Then, we note that, since transport is bal- 
listic, carriers can reach a given point inside the structure 
either directly from the contacts or indirectly after reflec- 
tion at the self-consistent potential barrier. Moreover, 
we can identify from which contact the carrier comes 
from. From these considerations it is easy to convince 
oneself that, in the low frequency limit of interest in the 
present paper, the distribution function solving the set 
Of Eqs. |)-(D is given by 

F(x,v x ,t) = 

' f (u, t) U < u < oo, v x > 0, < x < x m 

fh{u,t) U m < u < oo, v x < 0, < x < x m 

fo{u, t) U <u < U m , v x < 0, < x < x m 

fo(u, t) U m < u < oo, v x > 0, x m < x < L 

fh{u,t) U L <u < U m , v x > 0, x m < x < L 

, fh{u,t) U L < u < oo, v x < 0, x m < x < L 



(A2) 



where 



f {u,t) = F(o,+y/2 (u-U ),t) , (A3) 
f L (u, t)=F (L, -y/2(u- U L ),t) , (A4) 



with F (0, v x ,t) and F(L,v x ,t) being obtained through 
the boundary conditions in Eqs. (fJ) and (||). In the pre- 
vious expression U m = U m (t) and x rn = x m (t) refer to 
the value and location of the potential energy maximum, 
respectively, while Uq and Ul correspond to the values 
of the potential energy at the contacts (note that, ac- 
cording to the boundary conditions assumed here, these 
values do not flu ctuate). 

Since Eq. flA2| ) depends on the value of the maximum 
potential energy, to completely determine the distribu- 
tion function we need to derive the equation satisfied by 
the energy maximum. To this purpose, we first obtain 
an expression for the carrier density b y pe rforming the 
integral in Eq. (||) with the help of Eq. ( A2). After some 
algebra one can show that the carrier density is given by 



n(x, t) — 

n~ (U (x, t), U„ 
n+(U(x,t),U m 



; Lfo) h]) o < x < x m 

[fo, h]) x m < x < L 



(A5) 



where 



n-(U,U m -[f J L }) = 

" +0 ° du 



Um ^2{u - U) 



[fo(u,t) + f L (u,t)} 



du 



u y/2(u-U) 



2/oM) 



(A6) 



and 



n + (U,U m ;[f ,f L }) 



du 



Um y/2(u~U) 
u - du 



[fo(u,t) + f L (u,t)} 



u ^2{u-U) 



2/i(u,t) 



(A7) 



In the previous equations the square brackets on the l.h.s. 
mean a functional dependence. Now, by multiplying the 
Poisson equation, Eq.(Q), by dU(x, t)/dx and integrating 
with respect to a;, it can be shown that the electric field 
E(x, i) is given through 



E(x,t) = 

( E~ (U(x,t),U m ; [fo,f L ]) < x < x m 
\ E+ {U(x,t),U m] [f J L }) 

with 

E-(U,U m ;[f ,f L }) = 
du 



(A8) 



Vu-U- y/u - U m \ [fo(u,t) + f L (u,t)] 

1 

+ J duy/u-U2f (u,t) I 24 (A9) 



and 
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E+(U,U m ;[f ,f L }) 
du 



Vu-U-y/u- U m [f (u,t) + f L (u,t)} 

(A10) 



+ J duy / u-U2f L (u,t) \ 23 



Finally, from the definition dU(x,t)/dx = E(x,t) and 
its integration with respect to the space coordinate, we 
arrive at the following inverse equations for the potential 
energy profile 



dU 



Uo W{u,u m -[f ,f L ]) 



valid for < x < 

U + (x.t) 



X 1 1 



and 



dU 



Ul E+{U,U m ;[f ,f L }) 



~ x — L 



(All) 



(A12) 



valid for x rn < x < L. 

From the previous expressions the location of the max- 
imum ene rgy p otential can be obtained by evaluating ei- 



ther Eq. flAllD or Eq. (|Al|) at x 



dU 



u 



E-(U,U m ;[f J L ]) 
dU 



Ul E+(U,U m ;[f ,f L }) 



(A13) 
(A14) 



By eliminating x m from the two resulting equations we 
derive a closed equation for the value of the potential 
energy maximum in the form 



L = 



dU 



E-(U,U m ;[f ,f L ]) 
dU 

Uh E+(U,U m ;[f ,f L }) 



U 

u„ 



(A15) 



Equation ( A15 ) constitutes a closed equation to deter- 
mine the value of the potential energy maximum, U m . 
Notice that it depends solely on the boundary conditions 
for the distribution function and the sample length. 

Once the value of the potential energy maximum is 
know n, o ne can d eterm ine its location th roug h cither 
Eq. ( |A13| ) or Eq. ( |Al^ ). Then, from Eq. (|A|) one ob- 
tains the explicit expression of the distribution function. 
In a similar way the explicit spatial dependence of the 
potential energy can be de termi ned b y substituting the 
value of U m in Eqs. (All) and (A12), and that of the 
electric field and carrier density by substituting U m in 
Eq. (|Aq) and Eq. (A.5), respectively. In this way, a com- 



plete analytical solution of the model presented in Sec 
is obtained. 

In particular, the electrical current, defined as 



III 



I(t) = 



+ OC 



dv x v x F(x,v x ,t) 



(A16) 



can be shown to be given by 

/(*)= / du[f L (u,t)- f (u,t)] . (A17) 

J U m 

Now, we are in the position to derive the corresponding 
expressions for the transport and noise properties. 



1. Transport properties 

The average steady-state transport properties can be 
computed directly from the equations presented above by 
simply substituting into them 



X"fn j 

fo{u,t) -> f (u) = f c (u - Co) , 
h(u,t)^f L (u) = f c (u-Z L ) 



(A18) 



with / c (e) given through Eq . (^) . In this way we arrive 
at the equations used in Sec. IV A 



2. Low frequency current noise properties 

The fluctuating properties of any of the quantities of 
interest can be evaluated directly from the equations de- 
rived above by just performing the corresponding per- 
turbation around the steady state, and by taking into 
account that the only source of fluctuations in the sys- 
tem is located at the contacts, and represented by the 
fluctuating term in Eqs. ([?]). In this paper we are in- 
terested on the low frequency current fluctuations. To 
compute them we perturb Eq. (A17) around the steady 
state thus obtaining 



SI(t)= _ du[Sf L (u,t)-5f (u,t)] 
Ju m 

[fL(U m )~7o(U m )] SU m (t) . 



(A19) 



where 5U m (t) represents the fluctuations of the potential 
energy maximum. In the above expression we distinguish 
two contributions to the current fluctuation, one coming 
directly from the contacts and the other coming indi- 
rectly through the self-consistent potential fluctuations. 
To express the dependence of the se cond contribution on 
the noise sources, we perturb Eq. ( |A15| ). In performing 
such a perturbation it is convenient to shift all the energy 
integration variables by an amount equal to U m . After 
some algebra one then arrives at the following expression 
for the fluctuations of the maximum potential energy: 

5U m it) = i / + dvrf L (u)Sf L (u,t) + 
A Ju L 

Y r + co 

du^ (u)Sf Q (u,t) , (A20) 

Uo 
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with 



1 1 f + °° 
A ==--=- + / duj L (u)f FD (u - 

-C/L JC/r, 



+ 00 



+ / duj (u)f FD (u -£ ) 



where we have used that / (w) = — f fd{ u — Co) an d 
= — /fz>( u — Here, we have defined 

7L(u)=7 > (^(w-^m)+7i(^(F m -u) , (A22) 
7o(u)=7 > N^(^-^m)+7o < («)^(^ m -") , (A23) 
with 



7>(u) 



^2 (u - C7) - ^2 (u - C/ m ) 



S(C/,C/ m ;[/ ,/ L ]) 3 



y/2(u-U)-y/2 (u-U m ) 



dU 



E(U,U m ;[f J L ]) 3 



7o=(«) 



-2^2 (u-U)dU 
u L E+ {U,U m -[f J L ]f 

2^2 (u-U)dU 
u„ E-(U,U m ;[J J L }) 3 



(A24) 

(A25) 
(A26) 



By substituting Eq. ( |A20| ) in Eq. ( |A19| ) for the current 
fluctuations we finally obtain 



5I(t) 



dwy L (u)5f L (u,t) 



+ 00 



+ / dwfo(u)5f (u,t) 
Juo 



with 



7 L (m) = 9{U - U m ) + ^7i(lt) , 

7o(w) = —6{u - U m ) + ^7o(w) 



where 



Uc(U m -tL)-f c (U m -Zo)\ 

A 



(A27) 
(A28) 



(A29) 
(A30) 



(A31) 



Now we are in the position to compute the low frequency 
current spectral density, defined as 



5/(0) = / SI(0)SI(t)dt 



(A32) 



By substituting Eq. ( A27 ) into Eq. ( |A32j ) we arrive at 



5/(0) 



+oo 

Ul 
+oo 

Uo 



dwy L (u) 2 f FD (u -&) + 



du~fo(u) 2 f F D(u - Co) , 



(A33) 



(A21) 

where we have used that from Eq. (|^) one has 



2/ 5f a (u,t)Sf a ,(u',t>)dt' = 

J — oo 

fFD (u-USa,a>S(u-u')d(t-t') 



(A34) 



From Eqs. ( A24 )-( A32 ) it is str aightforward to arrive at 
the equations used in Sec. IV Bj . 
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